Width difference of the Bg mesons from lattice QCD 
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We extend our previous studies to calculate the B meson B parameters Bb and Bs on the lattice, and present an update of the 
results for the Bs meson width difference. We perform an extensive study of systematic errors in the quenched calculation of the 
B-parameters, and find that the systematic errors are in good control using the NRQCD action for heavy quark. We also report our 
preliminary results from unquenched simulations. 



1 Width difference AFs^ 

The width difference AF in the Bg — Bs system opens a 
new possibility to measure the CKM angles and to search 
for new physics, once it is measured and found to be siz- 
ableEl. Theoretical prediction of AF^ is, therefore, desir- 
able, and the most reliable calculation has been obtained 
using the Heavy Quark Expjansion under an assumption 
of the quark-hadron duality and a recent summary is 
given ina. In this calculation, nonperturbative inputs are 
necessary for the B meson B parameters Bb and Bs^ for 
which the lattice QCD may provide first-principles calcu- 
lation starting from the QCD lagrangian. In this report, 
we present an update of the lattice calculation of Bb and 
Bg by the JLiXID collaboration. Our previous works are 
published inii-Q. 

The main contribution to AF^^ comes from the cc 
final state, to which both B^ and B^ can decay, and the 
theoretical expression obtained using the Heavy Quark 
Expansion ° was given by Beneke et al. asHu 



X [cl{z){Ol)+cs{z){Os) 
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where Cl(z), cs{z) and Ci/miz) are known functions of 
z = rri^/ml. The four-quark operators Ol = h{l — 
Ibhfisbil - ^shfiS and Os = &(1 - J5)sb{l - 75)s are 
evaluated for Bs and Bg mesons as initial and final states 
respectively. Normalizing with the total width F^^ , one 
arrives at 



AF\ _ 16Tr^B{Bs Xev) PbMb, 
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" As pointed out by Falk at the conference, the operator product 
expansion could be dangerous for this quantity, because the energy 
release is not large enough for the cc final states. The quark-hadron 
duality could also be questionable, as the final states are saturated 
by a limited number of exclusive modes. Unfortunately, no rigorous 
method is known to estimate the violation of these assumptions. 
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where BB{rnb) and Bs{mb) are i?-parameters defined as 
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is a ratio of matrix elements of heavy-light axial current 
and pseudo-scalar density 0. |^ Other notations may be 
found in El. We note that for (AF/F)b^ the matrix ele- 
ment of Os gives dominant contribution while the effect 
of Ol is neghgible. The 1/m correction Si/^ is expressed 
in terms of matrix elements of higher dimensional Ai?=2 
operators, and is found to be sizable (and destructive) 
by an evaluation using the vacuum saturation approxi- 
mations. 

To obtain the prediction for (AF/F)b, , nonpertur- 
bative inputs are necessary for the decay constant /b^ 
and the S-parameters Bb and Bs- The quenched lattice 
calculation of the B meson decay constant is quite stable 
over years in the sense that results from many groups us- 
ing different lattice actions agree, and a recent summary 
is fss = 195(20) MeVcl. Several groups are recently per- 
forming unquenched calculations that include the effect 
of light quark loops. Although the calculation is compu- 
tationally much demanding, the unquenched QCD sim- 
ulation is being realistic on recent dedicated machines 
or on commercial supercomputers. Recent results sug- 
gest that the decay constants become higher with the 

The same notation Ti{mi,) is used for a different quantity in^. 
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Figure 1: Quenched results for Bsimi,) as a function of 1/Mp at 
/3=5.9. Results with four different methods are plotted. 



dynamical quark effect, and the most recent summajcy 
by Bernard at Lattice 2000 is /s, = 230(30) MeV 0, 
which we use in this work. 



2 Lattice calculation of B parameters 

Present work is ao-extension of previous calculations by 
Hiroshima arWPtfEl, and preliminary results are already 
reported in&O. 

We use the lattice NRQCD action for heavy quark 
including all 1/M corrections consistently. Therefore, no 
extrapolation in the heavy quark mass is necessary in 
contrast to the work by APE collaborations, in which 
they employ the relativistic action for relatively light 
heavy quark and extrapolate the results to mi,. 

We perform a study of systematic errors on quenched 
{Nf=Q) lattices. We calculate the _B-parameters in four 
different methods that have different systematic errors 
coming from higher orders in the 1/M and expansions, 
in order to see the associated systematic uncertainties. ^ 
In addition, the simulations are repeated at three lattice 
spacings to see the systematic error depending on lattice 
spacing a. 

Results for Bs{mi,) are plotted in Figure from 
which we can observe that the disagreement among the 
different methods becomes larger for larger 1/M. That 
behavior is expected, because the 1/AI expansion is trun- 
cated at first order in our calculation. Nevertheless the 
result at the Eg meson mass (1/Mp ~ 0.19 GeV~^) is 
under reasonable control. A similar plot is shown in Fig- 
ure ^ for three lattice spacings in the quenched approx- 
imation. We find a good agreement among the results 
from different lattices, which indicates that the system- 

For a detailed description of the four methods, seeB 
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Figure 2: Bs(mf,) as a function of 1/Mp at three lattice spacings. 
The method 2 is used. Unquenched result is also plotted (filled 
circles) . 



atic error associated with the lattice discretization is well 
controlled. 

At this conference, we presented a new result from 
unquenched simulations. As one of the major projects of 
the JLQCD collaboration, we are performing a two-flavor 
(dynamical u and d quarks) QCD simulation with a non- 
perturbatively 0(a)-improved quark at P—5.2, Csw=2.02 
on a 20^x48 lattice, which corresponds to the lattice 
spacing a ~ 0.1 fm. A preliminary report of this sim- 
ulation may be found inE3. On this lattice, we are cal- 
culating the B-parameters and a very preliminary result 
is plotted in Figure ||. Although the result is slightly 
lower than the quenched data, the difference is statisti- 
cally not significant, and our conclusion at this stage is 
that the effect of dynamical quarks is not sizable for the 
-B-parameters. 

Our (preliminary) results for _B-parameters in the 
quenched (iV/=0) and unquenched {Nf=2) calculations 
are 
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0.83(3) (8) [Nf = 2) 
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The systematic errors are evaluated using the variation 
among four different methods as a guide. Using (^ the 
result for the width difference is obtained as 



AT 



= 0.097t[!;Sj35 ± 0-025 ± 0.020 ± 0.016, (9) 



where errors are from residual scale /i dependence re- 
maining in the calculation of G{z) and Gs{z) in (^, Jb^, 
Bs, and the uncertainty in the 1/m corrections, respec- 
tively. The central value 0.097 is somewhat smaller than 
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our previous estimate 0.107(26)(14)(17)B due to a smaller 
input for /s^ , which was previously 245(30) MeV. The 
effect of unquenching is less significant (—6%). The ex- 
perimental results from LEE and SLD were summarized 
by Boix at this conferenceEl as (Ar/r)^^ = 0.16j:[J;^^. 

One of the major systematic errors in (|^) comes from 
the uncertainty in fs^ , which may be avoided by ciinsid- 
ering a ratio {AT/ AM)b,, once AMg is measuredu. We 
obtain 
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where errors come from /i, Bs/Bb, and the l/m correc- 
tions, respectively. 

3 Comparison with other approaches 

There are two other approaches to calculate the B meson 
U-parameters on the lattice. 

One is the HQET approach, in which the b quark is 
treated as a static color source, and it is a naive limit 
M ^ oo of the NRQCD action. Gimmfiz and Reyes 
calculated Bb and Bs in this method Eil^Ej, in which an 
unqucnched calculation is also performed (but with an 
unimproved action). 

The other is the relativistic approach, in which the 
relativistic lattice action is used for heavy quark. In or- 
der to avoid large discretization error growing as a power 
of aM, one cannot simulate the b quark and has to treat 
relatively light heavy quark with mass around 1-2 GeV. 
An extrapolation in the heavy quark mass is then at- 
tempted. The crucial point for this method is to find 
a heavy quark mass region where M is small enough so 
that the discretization error of 0((aA/)^) is under con- 
trol, and at the same time M is large enough to justify 
the use of heavy quark expansion in the extrapolation. In 
general it may not be easy to identify such a region with 
confidence, and the applicability of the method should 
be studied carefully depending on the quantities to cal- 
culate. Calculation of the iijparameters in_this method 
has been presented by APEB and UKQCdIIB collabora- 
tions. 

In Figure || we compare our results with other ap- 
proaches. While our data cover the physical region 1/Mp 
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Figure 3: Bs{mij) from different approaches. Our results are 
shown by open symbols {Nf=0) and solid circles {Nf=2). 



= 0.05-0.3 GeV-^ the HQET results are plotted at the 
static limit 1/Alp — 0, and the relativistic results are 
shown in the charm quark mass region l/Mp = 0.4- 
0.5 GeV~^. Although we did not attempt to fit the 
whole data with a single curve, it seems that the results 
of different approaches agree reasonably well. It suggests 
that the lattice calculation of Bs is in good shape, even 
though the individual method might have its own sys- 
tematic uncertainty. 

Besides the calculation of B-parameters, the APE 
collaboration proposed to express AF b^ using AMd as a 
normalizationt 
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and quoted a much smaller result 0.047 ± 0.015 ± 0.016. 
if is a known factor, and an experimental value may be 



used for (jb AmJ^^)- In the ratio ^ „ 

of systematic error should cancel in the lattice calcula- 
tion, and the result has much smaller uncertainty. The 
problem in this normalization is, however, that the large 
uncertainty enters implicitly through the (ratio of) CKM 
elements iVfs/Virfp, for which they use a value obtained 
from a global fit of the CKM elements. We do not take 
this strategy, because the several sources of (theoretical) 
systematic errors are hidden in the CKM fit, and it makes 
the estimate of systematic uncertainties less transparent. 
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4 Requirements for further improvement 

In order to improve the accuracy of the prediction (^), 
better determination of fs, is important. For the un- 
quenched lattice calculation, more work is necessary to 
reach the situation in the quenched case, where consis- 
tency is checked among several different approaches by 
many groups. The similar work is necessary for the B- 
parameters, for which only a few unquenched calculations 
are available including ours. 

In (H), it is evident that the 1/m correction gives 
a large negative contribution to the width difference. 
Therefore, its reliable evaluation is quite important. At 
present, it is estimated using the vacuum saturation ap- 
proximation, and thus nonperturbative calculation on 
the lattice is desirable. Although it is hard to deter- 
mine the (perturbative) matching coefhcients for higher 
dimensional operators, the lattice calculation should still 
be useful to have an idea how reliably the 1 /m corrections 
are estimated. 

5 Conclusions 

The JLQCD collaboration started a lattice calculation 
of Bb and Bs including the effect of sea quarks. A pre- 
liminary result shows that the quenching effect is not 
substantial for these quantities. 

We show that the calculations with several differ- 
ent methods on three lattice spacings show a reasonable 
agreement in the quenched approximation. It indicates 
that, using the lattice NRQCD, the systematic errors are 
under control. 

The results for (Ar/r)^^ still have large uncertainty. 
Better estimation of the nonperturbative inputs as well 
as the 1/m corrections will be necessary to improve the 
accuracy. 
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